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Using a fermionic renormalization group method we present a simple real 
space picture of the strong influence an impurity has on the electronic proper- 
ties of a Luttinger liquid. We compute the flow of the renormalized impurity 
potential for a single impurity over the entire energy range - from the mi- 
croscopic scale of a lattice- fermion model down to the low-energy limit. We 
confirm that low energy properties close to the impurity are as if the chain 
is cut in two pieces with open boundary conditions at the end points, hut 
show that this universal behavior is only reached for extremely large systems. 
The accuracy of the renormalization group scheme is demonstrated by a di- 
rect comparison with data obtained from the density-matrix renormalization 
group method. 

P ACS numbers: 71.10.Pm, 73.22.-f 73.63.-b 



1. INTRODUCTION 

The detailed theoretical understanding of the low-energy electronic prop- 
erties of one-dimensional interacting electron systems obtained within the 
last fifty years has mainly been reached by mapping microscopic models onto 
an effective field theory using the very powerful technique of bosonization.El 
Bosonization, renormalization group (RG) investigations,! and the use of 
conformal field theorytl have led to the concept of Luttinger liquid (LL) be- 
havior which unifies the low-energy physics of a wide class of models. In the 
mapping process terms which are expected to be irrelevant in the low-energy 
limit are neglected and the energy scales of the microscopic models get lost. 
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Therefore a direct investigation of lattice models is necessary to obtain in- 
formations about the energy scales at which the typical LL power-laws and 
thus universal behavior can be expected. This is of special importance in 
connection with the interpretation of experimental data. Such an inves- 
tigation might furthermore lead to a more intuitive understanding of the 
mechanisms involved in the physics of low-dimensional correlated electrons. 
Here we present such a direct fermionic discussion using analytical and nu- 
merical techniques focusing on the drastic influence a single impurity has on 
the low-energy physics of LL'sMlfflffi 

Within bosonization an impurity with scattering amplitudes V^k', (k 
denotes the wave vector) is modeled by forward and backward scattering 
only, i.e. k and k' are restricted to ±&f, where kp is the Fermi wave vector. 
Both scattering^ channels are furthermore decoupled. A perturbative bosonic 
RG calculation^ and a boundary conformal field theory analysis^ led to the 
following picture: For a chain of spinless fermionscl with repulsive interac- 
tions (LL parameter K p < 1) the backscattering amplitude Vb presents a 
relevant perturbation which grows as A^ p_1 when the flow parameter A is 
sent to zero. This leads to a breakdown of the perturbative analysis. The 
scaling behavior can be traced back to the power-law singularity of the 2fcp 
density response function in a LlJ^'S On the other hancl^ a weak hopping t w 
between the open ends of two semi-infinite chains is irrelevant and scales to 
zero as A k p ~ 1 . Assuming that the open chain presents the only stable fixed 
point it was argued that at low energy scales and even for a weak impurity 
physical observables behave as if the system was split in two semi-infinite 
chains with open boundary conditions at the end points^ Our main interest 
will be on the local spectral weight Pj{u) for lattice sites j close to the impu- 
rity and energies u close to the chemical potential p. For PjioS) a power-law 
suppression Pj(to) ~ \uj\ aB with the boundary exponent as = K^ 1 — 1 which 
only depends on the interaction strength and shape and filling of the band, 
but not on the impurity parameters, was predicted.^ Within the bosonic 
field theory the above assumption was verified by refermionization,! quan- 
tum Monte Carlo calculations ,H30 and the thermodynamic Bethe ansatz.H3 
Although this gives a consistent picture there are three important questions 
which have not been answered using field theory: 

• In a more elaborate RG procedure the growing backscattering am- 
plitude feeds back into the flow of all the other scattering channels 
Vk,k' neglected in the field theory. Does this lead to other stable fixed 
points? Does the above scenario hold for microscopic lattice models? 

• Provided the scenario holds what is the scale on which the universal 
open boundary fixed point (BFP) physics can be observed? 
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• Is it possible to obtain a simple physical picture of the "splitting" 
mechanism? 

To numerically investigate the first two questions exact diagonalization 
(ED) and the density-matrix renormalization group (DMRG) method were 
applied to the lattice model of spinless fermions with nearest neighbor in- 
teraction and a single impurity.l3'E3oJl3 The expected scaling. was confirmed 
for both weak impurities and weak hopping using ED data.0 However, due 
to the limited system size it was impossible to go beyond the perturbative 
(in either Vb or t w ) regime. Later it was claimed that the full flow from a 
weak impurity to the BFP was successfully demonstrated using DMRG,Ei'lll 
although this strong statement is not really supported by the numerical data 



presented. The smallest temperature discussed in Ref. 14 corresponds to a 
system of around 300 lattice sites and the largest system considered in Ref. 
|i~3| was N = 52, while in Ref. |i~5| it was shown that N ~ 10 2 lattice sites are 



clearly not enough to exclude an asymptotic behavior not governed by the 
BFP, even if one starts out with a fairly strong impurity. 

An attempt to answer the above questions requires a method which is 
non-perturbative in the impurity strength, does give the LL power-law di- 
vergence in the density response function, but is not limited to systems of a 
few hundred lattice sites. We here use a functional RG method which has 
recently been introduced as a new powerful tool in the theory of interacting 
Fermi systems.EjOta Details of the derivation of the flow equations are given 
and the RG scheme is applied to the spinless fermion model with site or hop- 
ping impurities. In the extended numerical analysis of the flow equations 
(see Sec. ^]) we mainly focus on the hopping immirity. This complements 
our earlier publication on the impurity problemJ13 In our RG approach the 
complete flow of the renormalized on-site energies and the renormalized hop- 
ping amplitudes from the microscopic energy scale down to the infrared fixed 
point is calculated. The flow equations are non-perturbative in the impurity 
strength while perturbative in the electron-electron interaction. We treat the 
full functional form of the renormalized impurity potential as generated by 
the flow, instead of replacing it approximately by the scattering amplitudes 
at the Fermi level. Computing the local spectral weight near the impurity 
we can confirm the flow to the BFP only if we start out with already fairly 
large impurities even for systems of up to 10 4 -10 5 lattice sites. This shows 
that for intermediate impurity and interaction parameters extremely large 
systems are required to reach the universal BFP. The quality of the approx- 
imations involved in our RG scheme is demonstrated by a direct comparison 
with essentially exact DMRG data for systems with up to N = 512 sites. 
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2. THE MODEL 

The one-dimensional lattice model of spinless fermions with nearest 
neighbor hopping amplitude t = 1, lattice constant a = 1, and nearest 
neighbor interaction U is given by 

H = - ( c iS+i + 4+i c i) + U H n i n i+^ (!) 

3 j 

in standard second-quantized notation. The boundary conditions we con- 
sider are discussed below. Here we mainly focus on the half filled band case 
for which the LL parameter 

(2) 

is analytically known from the Bethe ansatz£li As already indicated by this 
expression the model at half filling shows LL behavior only for —2<U<2. 
For U > 2 the groundstate displays charge density wave order and for U < 
—2 phase separation sets in. We restrict ourselves to repulsive interactions 
in the LL regime, i.e. to < U < 2. To the Hamiltonian H we either add a 
site impurity H s = Vrij or a hopping impurity = (1 — O( c i c jo+l +h.c). 
t' = corresponds to a vanishing hopping between the sites jo and jo + 1. 
Taking the thermodynamic limit or assuming open boundary conditions in 
Eq. (|T|) a small t' thus models the situation of a weak link between two 
decoupled chains. In the non-interacting limit the reflection coefficient for 
scattering of such impurities at fc_p = 7r/2 is given byll3 \R S \ 2 = V 2 /(A + V 2 ) 
for H s and |i4| 2 = (1 - t' 2 ) 2 /(l + t' 2 ) 2 for H h , which provides us with a 
measure for the strength of the bare impurity. 
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3. A HARTREE-FOCK STUDY 

As a first step it is instructive to consider the impurity problem within 
the Hartree-Fock approximation, before turning to the RG treatment. The 
impurity leads to Friedel oscillations in the non-interacting density nrofile 
{nj)o which for large \j — jo\ behaves as \R\ sin (2kp\j — io|)/|i — io[-" Sim- 
ilar oscillations are found in the matrix element (c^c J+1 )o. Thus both the 
Hartree potential = f/((nj_i)o + (nj+i)o) and the Fock "hopping cor- 
rection" U((cjCj +1 )o + c.c.) are oscillating and decay as l/\j — jo\- One then 
has to solve a (non-trivial) one-particle problem within such a non-local po- 
tential. This problem is interesting in itself - a special variant of it has been 
discussed in the mathematical literature£3 Taking into account the Hartree 
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term only, the resulting spectral weight for \u\ — *■ shows power-law behav- 
ior with an exponent which independently of the local part of the effective 
potential around jo is proportional to the amplitude U\R\ of the asymptotic 
oscillations £3 It can be shown analytically that this behavior is not changed 
when the Fock term is included. A mor e detailed discussion of this scat- 
tering problem will be given elsewhereB2l Thus solely due to the oscillations 
and the very slow l/\j —jo\ decay of the non-local effective potential already 
HF yields a power-law for the spectral weight, but with an exponent which 
not only depends on U, but via R also on the bare impurity strength. This 
contradiction to the bosonization prediction is not surprising as HF does 
certainly not contain the renormalization of the impurity strength. 

The extension of the HF study using self-consistent HF leads to un- 
physical results and thus cannot be used to gain further insight. The self- 
consistent iterative solution of the HF equations generates for all U > a 
groundstate which shows charge density wave orderta which is qualitatively 
incorrect since a single impurity cannot change bulk properties of the system. 



4. THE FERMIONIC RENORMALIZATION GROUP 

We now treat the problem using a fermionic functional RG approach. 
Various versions were applied recently to problems of strongly correlated 
two-dimensional electron systems. c3 Here we use the method proposed by 
Wetterichlll and Morrisjll where one introduces a cut-off parameter A in the 
free propagator G° cutting out degrees of freedom on energy scales less than 
A and derives an exact hierarchy of coupled differential flow equations for 
the one-particle irreducible vertex functions by differentiating with respect 
to A, where A flows from oo to 0. 

Before turning to the specific lattice model of spinless fermions with 
nearest neighbor interaction and hopping we will give a brief introduction 
to the method for a general lattice model of spinless fermions with inter- 
action and impurity free one-particle states \a) (e.g. Wannier states \j) 
or momentum states \k)), two-body interaction lEa^J^^ji^^S' 
and an impurity potential J2a,/3 Va,p c a c p- The exact flow equation for the 
selfenergy reads 



-i dG°> A {iui) 



dA 



1 - E A (^)G°- A (^) 



T a,i,P,s( iu; n , iu>i, w n , iui), 



(3) 



(5,7 
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where and S A are matrices with matrix elements G '^ and re- 
spectively, u n denotes a fermionic Matsubara frequency, and T is the tem- 
perature. The irreducible four point vertex T A p s obeys a flow equation 
with terms bilinear in T A and a term linear in the irreducible six point ver- 
tex. The initial condition is given by £^8°° = Va,/3- For spinless fermions 
the electron-electron interaction is renormalized only by a finite amount of 
order interaction squaredB Hence as our central approximation we replace 
the renormalized two-particle vertex to leading order in the interaction by 
the antisymmetrized bare coupling, i.e. r A 7(3(5 — * v a ^fi t $. With this ap- 
proximation S A becomes frequency independent. For the case of a model 
including a spin degree^ of freedom the flow of the two-particle vertex has to 
be taken into account.0 In order to simplify the remaining Matsubara sum 
in Eq. (||) we take the zero temperature limit and perform the RG flow with 
a frequency cut-off 

G°' A (iuj) = @(\uj\-A)G (iu J ). (4) 

The first factor on the right hand side (rhs) of Eq. (J3j) then involves products 
of a Dirac delta function and step functions. Using the relational 

5 e (lo - A)/ (0 e [u - A]) - 8{u - A) f 1 f(t)dt, (5) 

Jo 

where e is a broadening parameter tending to zero and / is a sufficiently 
smooth function, we obtain 

^ = ~ E E^^ 7 He M+ , (6) 

lu=±A 7,(5 



G A (iio) = \ G { \iuj) ' -£ A }> (7) 



-l 



-l 



where 



is the full propagator for the cut-off dependent selfenergy. The convergence 
factor e luj0+ with co = ±A is relevant only for determining the flow from 
A = oo down to some arbitrarily large finite Aq. For Aq much larger than 
the band width this high energy part of the flow can be easily computed 
analytically and yields the simple contribution j£ 7 \7,j3,7' One can then 
drop the convergence factor and continue the flow from Ao downwards with 
the new initial condition 



S i = ^ + lE^,A7' (8) 

7 
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Within this scheme ~ is the flowing effective impurity potential. In 
our RG procedure all the different impurity scattering channels correspond- 
ing to transitions from the one-particle states \a) to \(5) are coupled [see Eq. 
(|6|)] and we do not replace them by scattering amplitudes at the Fermi level. 
This has to be contrasted to the bosonization approach to the problem. The 
RG is a grand canonical method for which the chemical potential is fixed and 
the average particle number has to be determined from the Green function. 
To calculate observables, as e.g. the local spectral weight, we have to deter- 
mine the selfenergy at A = and invert the matrix Eq. (|7|) to obtain 
the Green function, i.e. solve the one-particle problem of a particle moving 
in the effective scattering potential E^°. This gives us a simple picture - in 
momentum or real space depending on the choice of one-particle states \a) 
- of the influence an impurity has on the electronic properties of a LL (see 
below). 

Because of the matrix inversion involved in calculating the rhs of Eq. 
(|6|) we were not able to analytically solve the flow equation for general pa- 
rameters. But it easy to show analytically that in the weak impurity limit 
Eq. (|6|) exhibits the scaling predicted from bosonization. For that purpose 
we work in momentum space and consider for k — k! different from a 
reciprocal lattice vector K. Then the term linear in E A presents the leading 
approximation in the expansion of G A on the rhs of Eq. (^) and we obtain 



fci,fc2 K 
1 _* 



J k2,k\ 



1 

« A - 61 



+ (A--A) 



(9) 



where = — jx with the one-particle dispersion e&, v(k) is the Fourier 
transform of the two-particle interaction, and iV is the number of lattice sites. 
We treat the Hartree-type and Fock-type terms separately. The Hartree-type 
term alone would lead to a selfenergy which depends on k — k! only. In the 
thermodynamic limit one obtains 



d v a 



H 



x£fc- fe , +A ,, ljfcl — 



*■ dki 

1 



1 



£A — t^k-k'+ki 
+ (A^-A) 



(10) 



For k = kp, k! = —kp the k\ -integration contains a singular contribution 
proportional to 1/A which leads to the power-law behavior with an exponent 
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proportional to v(2kp) as discussed below. The Fock-type term reads 



d_, 
dA' 



J k,k' 



1 r dko 



2tt j— tt 2tt 



v(k - k 2 ) 



xE 



1 



k2,k2+k'—k 



&A - ik 2 +k'-k 



«A - £k 2 

+ (A--A) 



Away from half filling the singular contribution ~ 1/A for k = kp, k' = —kp 
comes from k% ~ kp and is therefore proportional to v(0). If only the singular 
contributions are kept the differential equation for Ej^ _j tp reads 



dA S 



kp,—kf 



v(0)-v(2k F ) 



2nvp 



a Sa 



(12) 



where vf denotes the Fermi velocity. This leads to the scaling 



^k F ,-k F ~ 



I \ [v(p)-v(2kF)]/[2m> F ] 

A, 



(13) 



To leading order in the interaction the exponent [v(0) — v(2kp)]/[2irvp] is 
just 1 — K p which shows that our non-perturbative fermionic RG captures 



the power-law increase found in the perturbative bosonic RG. As Eq. (12) 
was derived by expanding the Green function G A in powers of the selfenergy, 
the scaling behavior Eq. (|l3|) can be trusted only as long as E A stays small. 
Eq. (^) also holds in the half filled band case since the two additional terms 
from the Hartree- and Fock-type contributions cancel each other. The scaling 
equation (|l3|) can similarly be derived for a continuum model. 



5. NUMERICAL SOLUTION OF THE RG EQUATIONS 

Numerically integrating the RG equations @ for finite systems of N 
lattice sites we can go beyond the perturbative regime. In each step of 
the integration we have to invert an N x N matrix. If we now specialize 
on the lattice model of spinless fermions with nearest neighbor interaction 
and hopping, treat the problem in real space, and assume open boundary 
conditions, then [G°(iu;)] 1 — E A is tridiagonal and Eq. (||) reads 

s=±l uj=±A 

= ^ E 4::lH' (15) 
oj=±A 
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Therefore the numerical effort is considerably reduced which allowed us to 
treat systems with up to 2 15 = 32768 lattice sites. Since the position of 
the impurity jo is chosen to be far away from the open boundaries they 
do not influence the behavior close to jo- We have confirmed this by also 
investigating systems with periodic boundary conditions, in which case we 
are limited to smaller systems of the order of 10 3 lattice sites. The initial 
conditions for the case of a site impurity are = U/2, Ejyjy = U/2, 

^fojo = V + U, and E • 9 = U for j ^ jo and not at the boundaries, while 
the other matrix elements are initially zero. We started the integration at 
Ao = 1000 convincing ourselves that further increasing Ao only leads to 
negligible changes in the results. For a hopping impurity we start with 
Eft = 17/2, = U/2, Ejo = u and s Ao , o+i = E Ao +i . o = ! _ t / 5 while 

the other matrix elements are again zero. 
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Fig. 1. SA=7° for a site impurity (circles and squares) with V = 1 and a 
hopping impurity (triangles) with t' = 0.5, both for U = 1. For a clear 
representation the nearest neighbor parts are multiplied by —3 respectively 
—2 as indicated in the legend. 



Fig. H shows typical results for and ^fj+i for U = 1, with a 

site impurity V = 1.5 in one case and a hopping impurity t' = 0.5 in the 
other, where both would lead to a transmission probability of |T| 2 = 0.64 
in a non-interacting system. Since S A is symmetric around jo and the ends 
of the chain are dominated by the open boundaries mainly the region 1 <C 
j < jo is shown. For symmetry reasons in the first case we take N = 1025, 
jo = 513 and the average number of electrons at fj, = U is (Np) = 512 
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so that we are slightly off half filling. In the latter N = 1024, jo = 512 
and (Np) = 512. At half filling a hopping impurity leads to the shown 
long range oscillatory behavior in but not in similar to the HF 

approximation. Furthermore both the diagonal and the nearest neighbor 
part show global shifts. For this shift by U is due to the initial condition 
and is exactly canceled by the chemical potential which at half filling is given 
by /j, = U. A change in is present already in HF and corresponds to 

an interaction induced broadening of the band. 
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Fig. 2. - Q as a function of A for a site impurity (same parameters as in 
Fig. 1) and different N. 
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Fig. 3. The effective hopping 1 — - 0+1 between sites jo and jo + 1 as 
a function of A for a hopping impurity (same parameters as in Fig. 1) and 
different N. 
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In Fig. |H we show - o for the site impurity as a function of A for 
different N on a log- linear scale. For finite N the flow is effectively cut 
off on a scale of the order of 1/N and to extrapolate to the thermodynamic 
limit we have to consider a sequence of N values. Obviously the renormalized 
potential at the impurity site remains finite in the limit N — > oo and the 
expected "cutting" of the chain does certainly not occur because a single 
on-site energy diverges, as one might guess if the bosonic RG is taken too 
literally. Singular behavior is only found in for momenta with k — k'~ 
±2/cp, which is associated with the long range oscillations in real space. 

Fig. H shows the effective hopping 1 — - o+1 for the hopping impurity 
between sites jo and jo + 1 as a function of A for different N on a log-linear 
scale. The effective hopping does not tend to zero for N — > oo as one could 
expect from a simplistic interpretation of the bosonization result. As we do 
not perform a rescaling in our RG analysis the weak hopping has to stay 
finite because we otherwise would suppress the transmission through the 
weak link on all energy scales. The hopping between jo and jo + 1 even 
increases upon integrating out degrees of freedom. Again the long range 
oscillations of the non-local effective impurity potential and not the scaling 
of a single on-site energy or hopping matrix element is the reason for the 
peculiar behavior of physical observables, as for example pj(w), discussed 
next. 
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CO 

Fig. 4. pj -i(uj) as a function of cu for a site impurity with the same para- 
meters as in Fig. 1. 



The local spectral function at the site jo — 1, pj —i(u), is presented in 
Fig. H for a site impurity with the same parameters as in Fig. ffl. The data 
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show a suppression of the weight for \u\ — > 0, i.e. for energies close to the 
chemical potential, as expected from bosonization. A similar suppression is 
found in the local spectral weight close to a hopping impurity. Each spike 
represents a 5-peak of the finite system. For large systems the suppression 
follows a power-law which similar to the HF description is related to the os- 
cillatory long range nonlocal effective potential. But while the HF exponent 
depends on the bare impurity strength we expect our RG method to lead to a 
power law decay described by the boundary exponent a^ G (U) independent of 
the the impurity strength, where a>g G (U) can be determined from the same 
RG analysis for an impurity free system by calculating Pj{uS) for j close to 
one of the boundaries. The occurrence of ctg G (U) can be shown analytically 
from Eqs. ([H] ) and (^) only for \V\ S> 1 or \t'\ -C 1. We shortly sketch the 
argument for the weak hopping case. If one introduces operators L and R 
which project on the one-particle states to the left and right of the weak link 
the resolvents G A (z) = (z — in Eqs. (14) and ( |i~5| ) can be expanded 

in the effective hopping matrix element /i^ 0J0+1 = — 1 + ^f j t 
anticipating that this quantity is small for all values of A. 



yA 



G LL (z) = LG A (z)L = [zL - Lh K L}- 1 + O 



yA 



One obtains 



and a correspond- 



ing result for RG A (z)R. If we denote the leading term as 



Gi L (z) n (and 



G K RR {z) 



and also expand LG A (z)R = [G ll(z)] Lh A RG rr(z) the equa- 



tion for ^f j +i reads 



dA j0 ' J °- 



1 ~ ^ ? A ^°' i0 

\ uj=±A 



/~iA 



iH 



jo, jo- 



(16) 



The correction terms to the rhs are of order (^f j +ij ■ We also allowed 
the Coulomb integral U' across the weak link to be different from U. The 
equations for all other nonvanishing matrix elements of S A , apart from higher 
order corrections, are equations for a chain split between the sites jo and 
jo + 1. Therefore it is obvious that local spectral functions Pj(co) for j in 
the neighborhood of jo show a power law suppression with the boundary 
exponent ag G (U). Once the equations for the split chain are solved the 
result for follows from Eq. (|i~6| ) by a simple integration. In the 

special case U' = the effective hopping is given by t' for all values of A. 
The case of a strong site impurity | V| 3> 1 can be mapped to a weak hopping 
problem with t' ~ 1/V hy projecting out the impurity site. 

It is advantageousEj to analyze the finite size scaling of the spectral 
weight W(N) at [i and jo — 1 instead of trying to fit a power-law to the 
fixed N data for pj _\((J) at small uj. If pj _i(u;) follows a power-law as a 
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function of frequency, we obtain a power-law with the same exponent in the 
N dependence of W(N). In the next section we investigate the behavior of 
the exponent as a function of system size and compare the data obtained 
within the RG scheme with essentially exact data from DMRG focusing on 
hopping impurities. 



6. LARGE N AND COMPARISON WITH DMRG RESULTS 

We have calculated the spectral weight W(N) at \i and on the site next 
to the modulated bond jo for different t', U, and N = 2 n . Applying the 
RG scheme we were able to go up to n = 15, while using the DMRG we 
were limited to n = 9. Since we are interested in the exponent of a possible 
power-law behavior we determine 

/at nn\ In [W (2-+ 1 )] - In [W (2^)] 

ai{N " 2 } " ln[2»+i]-ln[2«-i] ' 

If W(N) decays for ff^ooasa power-law, aj(N) converges to the respec- 
tive exponent. 
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Fig. 5. ai(N) as a function of N for £7 = 0.5 and different f': f = 0.8 
(circles), t' = 0.5 (squares), t' = 0.2 (diamonds), and t' = (triangles). The 
filled symbols are DMRG data, while the open ones have been obtained from 
the RG. The dashed-dotted line gives the exact boundary exponent . 



Fig. m shows oli{N) for U = 0.5 and t' = 0.8 (corresponding to |T| 2 = 
0.95 at U = 0, i.e. a very weak impurity), i' = 0.5 ([T| 2 = 0.64 at U = 0, 
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i.e. an intermediate impurity), and t' = 0.2 (|T| 2 = 0.15 at U = 0, i.e. a 
weak hopping) . For comparison we also calculated a b (N) for the lattice site 
next to an open boundary (t' = 0). The DMRG and RG data are parallel 
to each other, which in addition to the analytical arguments is a strong 
indication that our fermionic RG captures the essential physics. For t' = 
both methods produce the expected power-law behavior with boundary 
exponents ag MRG and a| G . ag MRG (A = 512) agrees up to 1% with the 
exact exponent a|f f» 0.1609 obtained using Eq. (Q) and af = K~ l — l£9 
As our RG is only correct to leading order in U, a RG (A = 16384) which 
effectively is equal to a RG (A = oo) deviates by roughly 6% from a e g. From 
the RG curves for small to intermediate bare hopping t' between jo and 
jo + 1 one can infer that a RG (A) will approach the impurity independent 
boundary exponent a RG for A — > oo, as predicted by bosonization. We 
have also determined af G (N) for a site impurity and three different values 
of V leading to the same U = transmission probabilities |T| 2 as for the 
impurity parameters discussed here (see above). A comparison shows that 
for the hopping impurity and a fixed |T| 2 the convergence to a B G is slightly 
slower. Already for the site impurity extremely large = 10 4 -10 5 were 
needed toexclude non-universal (impurity dependent) fixed points with some 
certaintyEj and here even longer chains are required. 




Fig. 6. aj(N) as a function of N for if = 0.5 (dashed lines) and as (A) for 
t' = (solid lines) for different U: U = 0.5 (circles), U = 1 (squares), and 
U = 1.8 (diamonds). Filled symbols are DMRG data, open ones RG results. 
The dashed-dotted lines give the exact U dependent boundary exponents 
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In Fig. [6| RG and DMRG data are presented for an intermediate impu- 
rity strength tf = 0.5 and t' = for different values of U. Due to higher order 
corrections in U, the difference between the RG and DMRG data increases 
with increasing U. For the case of a site impurity the speed of convergence 
of af G (N) to q§ g was increased considerably by increasing U (see Fig. 4 of 
Ref . |l~9| ) . In Fig. |6| we also find a tendency towards a faster convergence to 
a§ G if U is increased for a fixed if , but the situation is less clear since the 
af~ G (N) curves at large N cross a^ G (N) showing non-monotonic behavior 
and approach a§ G from above. 

Compared to the case of a site impurity we here find somewhat less clear 
indications that the boundary exponent and thus the BFP is approached 
independently of the bare impurity strength. Nonetheless the presented 
RG data lead us to conclude that the universality of the BFP holds also 
for a hopping impurity. Our data show that for a hopping impurity and 
especially for intermediate t' and U, which are experimentally most relevant, 
exceptionally large systems are needed to observe the universal BFP physics. 
For chains which are not long enough a strong system size dependence of 
experimentally extracted exponents must be expected. For both types of 
impurities it is obvious that a DMRG study alone, which is limited to a 
few hundred lattice sites, cannot give a definite answer to the question if 
the BFP is indeed universal.EJ'Ej To judge the quality of the approximations 
made to set up our RG scheme the DMRG data are on the other hand very 
useful. 



7. SUMMARY 

In summary, we have applied a functional RG method to a lattice 
fermion model in order to analyze how and on what scale a single impurity 
in a Luttinger liquid "cuts" the chain. The method is perturbative in the 
electron-electron interaction but non-perturbative in the impurity strength. 
The results have been checked against essentially exact DMRG data for fi- 
nite systems, with good quantitative agreement for small U. Computing the 
local spectral weight near a hopping impurity we found a clear indication 
of the expected BFP physics only for intermediate to strong hopping impu- 
rities even for systems of up to N = 2 15 latticesites. Having in mind the 
results of our previous study for site impuritiestJ we nevertheless conclude 
that universal power-law behavior with a boundary exponent depending on 
the electron-electron interaction and not on the impurity strength also holds 
for hopping impurities. The functional RG method, where the entire func- 
tional form of the impurity potential is renormalized, provides a clear real 
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space picture of the cutting mechanism for both site and hopping impurities: 
while the local impurity strength is not renormalized much, an oscillating, 
slowly decaying non-local potential develops which at low energies acts as 
an effective barrier. 

It turned out that extremely large systems are required to reach the 
universal regime, even for intermediate impurity and interaction strengths. 
Hence only special mesoscopic systems, such as very long carbon nanotubes, 
under very favorable experimental conditions are suitable for experimentally 
observing the impurity induced universal boundary physics. A trend towards 
cut chain behavior is certainly easier to observe. In Fig. || e.g. a strong 
suppression of the local spectral weight can be seen already for a chain with 
N = 1025 lattice sites, a system size for which aj(N) is still far away from 
the universal boundary exponent as- Functional RG methods applied to 
realistic models should be useful to compute the non-universal behavior at 
intermediate scales, which is more easily accessible to experimentalists. 



ACKNOWLEDGMENTS 

We would like to thank W. Apel, P. Durganandini, P. Kopietz, N. Shan- 
non, C. Wetterich, and especially M. Salmhofer and H. Schoeller for valuable 
discussions. U.S. is grateful to the Deutsche Forschungsgemeinschaft for sup- 
port from the Gerhard-Hess-Preis. 



REFERENCES 

1. For a review see: J. Voit, Rep. Prog. Phys. 58, 977 (1995). 

2. J. Solyom, Adv. Phys. 28, 201 (1979). 

3. A. Luther and I. Peschel, Phys. Rev. B 9, 2911 (1974). 

4. D. C. Mattis, J. Math. Phys. 15, 609 (1974). 

5. W. Apel and T. M. Rice, Phys. Rev. B 26, 7063 (1982); T. Giamarchi and H. J. 
Schulz, ibid. 37, 325 (1988). 

6. C. L. Kane and M. P. A. Fisher, Phys. Rev. B 46, 15233 (1992). 

7. S. Eggert and I. Affleck, Phys. Rev. B 46, 10866 (1992). 

8. D. Yue, L. I. Glazman, and K. A. Matveev, Phys. Rev. B 49, 1966 (1994). 

9. Here we will mainly discuss spinless fermions. 

10. K. Moon, H. Yi, C. L. Kane, S. M. Girvin, and M. P. A. Fisher, Phys. Rev. 
Lett. 71, 4381 (1993). 

11. R. Egger and H. Grabert, Phys. Rev. Lett. 75, 3505 (1995). 

12. P. Fendley, A. W. W. Ludwig, and H. Saleur Phys. Rev. Lett. 74, 3005 (1995). 

13. S. Qin, M. Fabrizio, and L. Yu, Phys. Rev. B 54, R9643 (1996). 

14. S. Rommer and S. Eggert, Phys. Rev. B 62, 4370 (2000). 

15. V. Meden, P. Schmitteckert, and N. Shannon, Phys. Rev. B 57, 8878 (1998). 



A Single Impurity in a Luttinger Liquid 



16. C. Wetterich, Phys. Lett. B 301, 90 (1993). 

17. T. R. Morris, Int. J. Mod. Phys. A 9, 2411 (1994). 

18. M. Salmhofer, Renormalization (Springer, Berlin, 1998). 

19. V. Meden, W. Metzner, U. Schollwock, and K. Schonhammcr, cond- 



mat/0104336| (2001). 

20. F. Woynarovich and H. P. Eckle, J. Phys. A 20, L97 (1987); C. J. Hamer, G. 
R. W. Quispel, and M. T. Batchclor, ibid. 20, 5677 (1987). 

21. We have neglected a phase shift which is irrelevant for the present discussion. 

22. M. Klaus, J. Math. Phys. 32, 163 (1991). 

23. K. Schonhammer and V. Meden, to be published. 

24. A. Cohen, K. Richter, and R. Berkovits, Phys. Rev. B 57, 6223 (1998). 

25. D. Zanchi and H. J. Schulz, Phys. Rev. B 61, 13609 (2000); C. J. Halboth and 
W. Metzner, ibid. 61, 7364 (2000); C. Honerkamp, M. Salmhofer, N. Furukawa, 
and T. M. Rice, ibid. 63, 5109 (2001). 

26. For N = 2 9 , a?^ RG (iV) was determined by fitting In [W (2™)] as a function of 
In [2™] by a spline and then taking the centered difference according to Eq. (fli 



